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ABSTRACT 

Results of recent observations of the Galactic bulge demand that we discard a simple picture of 
its formation, suggesting the presence of two stellar populations represented by two peaks of stellar 
metallicity distribution (MDF) in the bulge. To assess this issue, we construct Galactic chemical 
evolution models that have been updated in two respects: First, the delay time distribution (DTD) 
of type la supernovae (SNe la) recently revealed by extensive SN la surveys is incorporated into the 
models. Second, the nucleosynthesis clock, the s-processing in asymptotic giant branch (AGB) stars, 
is carefully considered in this study. This novel model first shows that the Galaxy feature tagged 
by the key elements, Mg, Fe, Ba for the bulge as well as thin and thick disks is compatible with a 
short-delay SN la. We present a successful modeling of a two-component bulge including the MDF 
and the evolutions of [Mg/Fe] and [Ba/Mg], and reveal its origin as follows. A metal-poor component 
(<[Fe/H]>~-0.5) is formed with a relatively short timescale of ~1 Gyr. These properties are identical 
to the thick disk's characteristics in the solar vicinity. Subsequently from its remaining gas mixed with 
a gas flow from the disk outside the bulge, a metal-rich component (<[Fe/H]>~+0.3) is formed with 
a longer timescale (~ 4 Gyr) together with a top-heavy initial mass function that might be identified 
with the thin disk component within the bulge. 

Subject headings: Galaxy: bulge — Galaxy: abundances — Galaxy: evolution — stars: abundances 



1. INTRODUCTION 

The Galactic bulge is classified as a boxy bulge that is 
associated with bars and is l ikely generated through disk 
instability processes (e.g ., iKuiiken fe Merrifieldl Il995t 
I Bureau fe Freeman! Il999t ). Its kinematics and surface 
brigh tness profile are indeed shown to be cy lindrical rota- 
tion (|Howard et alJl2009b IShen et al.ll2O10l) and be near- 
exponential rather than r 1 / 4 (|Kent et al.lll991l ). respec- 
tively as evidence for boxy bulge. This view suggests that 
the formation of the Galactic bulge is closely connected 
to disk evolution, resulting in a longer timescale for its 
formation than that expected for classical bulges, con- 
sidered to be merger products as a result of hierarchical 
galax y formation in the cold dark matter (CDM) Un i- 
verse (jAguerri et al.ll20"01t IScannapieco fe Tisseral [2003). 
Thus, this views favors a recent finding of a large age span 
of ~2-13 Gyr amon g microlensed turn-off bulge stars 
(jBensbv et al.l [201lT ). though the color-magnitude dia- 
grams of the Galactic bulge show no clear si gnature of the 
presence of a young stellar population (e.s. IZoccali et all 
[20031: iClarkson et alJ[20M) . 

In general, it is expected that a disk instability forming 
the bulge induces a vertical mixing, which leads to eras- 
ing a metallicity gradient along a minor axis. This con- 
trad icts the observed re sult showing a clear [Fe/H] gradi- 
ent (jZoccali et alj |2008). This problematic issue is solved 
by introducing two components, i .e., thick and thin 
disks, as the origin of Galactic bulge (|Bekki fe Tsujimotol 
I2011al) . Their scenario is as follows: The first disk is 
disturbed by an ancient minor merger, which induces a 
vertical growth of the disk and transforms it into a thick 
disk, and subsequently the thin disk starts to form with 
an accompanying bar formation in the central region. In 
such a two-component disk, a vertical mixing induced by 



a bar buckling functions incompletely in a sense that the 
high latitude region in the thick disk is not well mixed. 
In the end, the resultant bulge shows a vertical metal- 
licity gradient as well as a cylindrical rotation, both of 
which are compatible with the observed results. 

Two co mponents of the Galact ic bulge were first pro- 
posed by iBabusiaux et "ail (|2010f ) in terms of chemistry 
and kinematics. They find the metal-poor population 
kinematically associated with an old spheroidal or a 
thick disk and the metal-rich population with a bar- 
like kinematics. Furthermore, very recently, a double- 
peaked metallicity distribution (MDF) of bulge stars 
has been repo rted with two d ifferent tracers, i.e., red 
clump giants (|Hill et al.l |2011[) and microlensed dwarf 
stars (Be nsbv et al.l I2011T) . Two results basically point 
to the same conclusion that the metallicity of one peak 
is metal-poor at [Fe/H]~ -0.6 - -0.3 and the other metal- 
rich at [Fe/H]~+0.3. 

In view of chemical evolution, the MDF should be 
examined together with the time evolution of the ele- 
mental abundance pattern generally described as [X/Fe] 
vs. [Fe/H]. One aspect of the converged knowledge on 
this for the bulge is that the [a/Fe] ratio decreases 
with an increasing [Fe/H ] for the m e tal-ri c h regime (e.g., 
Alves-Brito et alj 12010: Iffill et al.l 120111: iBensbv et all 
20111 ). which bears witness to the presence of contri- 
bution from type la supernovae (SNe la) to chemi- 
cal enrichment in the bulge. Here we should high- 
light recent results regarding the delay time distribu- 
tion (DTD) of SNe la yielded by the st udies on the SN 
la rate in distant and nearby galaxies jM annucc i et al. 
20061: ISulliyan et alj|2006t iTotani et~aT]|2008l:lMaoz et al. 
2010. 20111 ). Their findings dramatically shorten the SN 
la's delay time, compared with its conventional timescale 
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of ^1 Gyr, which will have a significant impact on Galac- 
tic chemical evolution. 

This renewed picture of a SNIa clocking should be tied 
up with another nucleosynthesis clock, the s-process op- 
erating in an asymptotic giant branch (AGB) star. These 
two different delayed-timings to release the nucleosyn- 
thesis product will create variation in stellar abundances 
among the elements such as Mg, Fe, and Ba, from which 
we can decipher the evolutionary pathway to the present 
through detailed modeling. In this paper, incorporating 
the new SN la DTD and AGB yield into the Galactic 
chemical evolution (GCE) model, we explore the chem- 
ical evolution of the Galactic bulge, highlighting its two 
components. To strengthen the rigidness of the model, 
we first complete the reproduction of the chemical fea- 
ture of the Galactic thin and thick disks. 

2. TWO NUCLEOSYNTHESIS CLOCKS 

Chemical enrichment at an early stage is promoted by 
the product from short-lived massive stars through type 
II SNe (SNe II), and subsequently with a time delay, 
heavy elements from SNe la and AGB stars come to an 
entry. As a result, the relative contribution from SNe 
la or AGB stars to the ISM in comparison with SNe 
II, which is imprinted in individual stellar abundances 
of long-lived stars, can be utilized as an age-dating for 
the stars. The updates of the input for these delayed 
contributors are briefly described. 

2.1. SNe la 

It is difficult to estimate the evolutionary timescale 
of SN la progenitors from purely theoretical arguments 
due to difficulties in identifying the binary companions 
of SN la's progenitors. Therefore, the break in the 
[a/Fe] ratio seen in the solar neighborhood stars due to 
the switchover of the major Fe source from SNe II to 
SNe la at [Fe/H] ~-l had long been the only informa- 
tion available for deducing the delay time of SNe la by 
theoretically estimating the elapsed time until [Fe/H]~- 
1. This assessment indeed leads t o a considerably long 
delay time of ~1-1.5 Gyr (e.g., iPagel fe Tautvaisienel 
Il995t lYoshii et all 1 19961 ). However, recently we have en- 
tered a new phase in our understanding of SN la de- 
lay time. Its distribution, i.e., the DTD, is revealed by 
the survey of SN la rate in extragalaxies, the individ- 
ual galaxy ages of which are deduced. Surprisingly, it 
turns out that young progenitors for SNe la are dom- 
inant l|Mannucci et all 120061: [Sullivan et all |2006|L and 
that the DTD is proportional to tj 1 for the approxi- 
mate nmgejafjll Glyr < ideiay <10 Gyr (|Totani et al.l 
l2008tlMaoz et alJl2010t l. This implies that about 70% of 
SNe la explodes with a time delay within 1 Gyr. Such a 
form of DTD can be predicted by the theoretical mod- 
els for the progenitor of SNe la: a double-degenerate 
(DD) scenario, i.e., merging of double C+O white dwarfs 
with a combine d mass surpassing the Chandrasekhar 
mass limit (e.g- lYungelson fe Lividl2000l : lGreggidT 2005) 
or a single-degenerate (SD) scenario, i.e., accretion of 
hydrogen-ri ch matter via mass transfer from a binary 
companion (jHachisu et al.l l2008) . In this study, we apply 
this new SN la DTD to Galactic chemical evolution, and 
for its normalization we assume the mass of the primary 
component of binary, which eventually produce SNe la 
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Fig. 1 . — Comparison of the theoretical [Ba/Fe] ratio expected for 
AGB stars as a function of [Fe/H] with the observed abundances of 
Post-AGB stars (red circles), Ba stars (blue squares), and Carbon 
stars (crosses). Solid lines denote the ratio of Ba yield synthesized 
in AGB stars to stellar metallicity used for model calculations. The 
Ba yield corresponding to dashed line showing a better fit to the 
observation is also applied to calculations (see Fig. 3). Two dotted 
curves represent the maximal and minimal pre dictions from AGB 
models with a range of 13 C-pocket efficiencies (Busso etaDHSS3). 



through accretion or merging, in the range of 3 — 8 Mq . 
In the end, the SN la rate is calculated so that the frac- 
tion of the stars that end up with SNe la for this mass 
range is 0.08, which is distributed with a DTD oc tj elay 
for 0.1 < tdeiay < 10 Gyr. Its fraction has been obtained 
through previous works as well as recheck by this study. 
In §3, we will see a fatal problem that we have noticed 
in a previous approach to try to determine the SN la 
DTD within the Galaxy thanks to advancement of the 
observed data. 

2.2. s -processing 

Low-mass AGB stars (M <3M Q ) release the s-process 
elements during th e thermally pulsing AGB phase (e.g., 
iGallino et a.1.1 fl998h . Mainly due to large uncertainties 
in convective mixing and 13 C-pocket efficiencies, the s- 
process nucleosynthesis allows a wide range for the level 
of a possible production. On the other hand, abun- 
dances of the surface of AGB stars can be directly com- 
pared with the nucleosynthesis results. Here we focus 
on the clement, Ba. Figure 1 shows the observed abun- 
dances of Post-AGB stars, Ba stars, and Carbon stars 
l|Kappeler et aLlfeOlll references therein). Note that all 
observed data resid e in the theoretically allowable range 
(|Busso et al.1 l2001h denoted by dashed lines. Here we 
determine the best empirical Ba yield as a function of 
stellar [Fe/H] so as to reproduce the Ba evolution of thin 
disk stars within an observationally acceptable range in 
terms of AGB abundances. For [Fe/H]>-1, the Ba yield 
given by dashed line which is fully consistent with the 
data fails to reproduce the chemical evolution of the thin 
disk. Alternatively, we adopt the Ba yield for stars with 
a mass of 1.5-3 Mq corresponding to the solid line which 
is rather close to the upper envelope of the observed 
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Fig. 2. — Predicted features of chemical evolution in the 
thick disk, compared with the observed quantities. Top panel: 
The [Mg/Fe1- [ Fe/Hl diagram. The obser ved d ata are taken from 
IBensbv et all (120051) . iRuchti et al.l i20Wi , and I Venn eTaTI (12001 1. 
denoted by filled circles, small circles, and small crosses, respec- 
tively. Mid dle panel: The [Ba /Mg]- [Fe/H] diagram. The [Ba/Mg] 
ratios from Venn et al. ( 2004) are shifted by +0.15 dex to broadl y 
adjust their average to the level of data from Bcnsby et al. (2005). 



Fig. 3. — Predicted features of chemical evolution in the thin 
disk. In the models, chemical evolution is assumed to start from 
the remaining gas after the formation of the thick disk, and the 
results are shown by solid curves. The dashed curve is the result 
calculated with the Ba yield, the metallicity dependence of which 
corresponds to th e dashed line in Fig.l. The [Ba/Mg] ratios from 
IVenn et all (120011 are shifted by +0.15 dex as in Fig. 2. 



[Ba/Fe]-[Fe/H] plane. 

3. CHEMICAL EVOLUTION OF DISKS 

Solar neighborhood stars are a mixture of stars be- 
longing to different Galactic components, i.e., halo, thick 
disk, and thin disk. A growing understanding of the 
kinematic properties of individual stars thanks to Hip- 
parcos data has enab led us to precisely se lect thick and 
thin disk stars (e.g., iFeltzing et al.l [20031) . It then re- 
veals that the observed break in [a/Fe] is seen among 
thick disk stars (Fig. 2), and that all thin disk stars fol- 
low a decreasing [a/Fe] trend (Fig. 3). Accordingly, our 
previous view that some fraction of thin disk stars pop- 
ulates a plateau of [a/Fe] with respect to [Fe/H] has to 
be discarded. In this study, we examine the evolution 
of [Mg/Fe] and [Ba/Mg] (not [Ba/Fe] so as to prevent 
the effect of s-processing from being hidden by SN la 
contamination). In addition, we focus on the chemical 
evolution for [Fe/H] since the origin of metal-rich 
disk stars should be assessed with an extra evolution fac- 
tor such as either stellar migration (iRoskar et al.l 120081: 
Bekki fc Tsuiimo"tol 1201 lbl) or bulge winds |Tsuumotol 
2007t iTsuiimoto et al.ll2010l) . 

We should note that Ba is synthesized through not 
only s-process but also r-process. The theoretical inter- 
pretation of abundance data on very metal-poor stars 
implies that the mass range for the r-process is 8-10 Mq 
([Mathews et al.lll992Ulshimaru fc Wanaidll999D as iden- 
tified with the col lapsing O-Ne-Mg core (Wheele r et al.l 
11998ft. 20-25 MrT, (ITsuiimoto et al.l l2"000l). or 12-30 M 
( Cescutti et al.ll2006ft with heavy weighing on lower mass 



progenitors. Here, as the site for r-process, we adopt the 
mass range of 2 -25 and the Ba yield deduced by 
ITsuiimoto et al.l ((2000). We do not have to be exact in 
the mass range because its variance influences only the 
behavior of r-process enrichment at a very early epoch, 
outside the focus of this study. 

Since the present GCE models are essentially the 
same as those adopt e d in our previous paper (e.g., 
iTsujimoto etaLl 120101 : [Bekki fc Tsuumotd l2011bh . we 
briefly describe the models. The basic picture is that 
the Galactic thin and thick disks were formed by the gas 
infall from outside the disk region. For the infall rate, we 
apply an exponential form with a timescale Tj n . Taking 
into account the relatively rapid formation of the thick 
disk and the presence of a G-dwarf problem for the thin 
disk, we assume a rather short timescale of Tj n =0.5 Gyr 
for the thick disk, and a much longer timescale of Tj n =5 
Gyr for the thin disk. The metallicit y Z- ln of infalling gas 
is assu med to be pre-enriched (see iBekki fc Tsuiimotol 
I2011bft . For the thick disk, we set [Fe/H] = -1.3 since 
the MDF of thick disk sta rs shows a sharp increas e from 
[Fe/H]— 1.3 to the peak dWvse fc Gilmord Il995l) . For 
the thin disk, the metallicity is determined by an impli- 
cation from the cos mic evolution of damped Lya systems 
(jWolfe et al.l 12005). the metallicity of which is around 
[Fe/H] = -1.5 at the epoch of thin disk formation. Thus, 
in our models an infalling gas for the thick disk is some- 
what more enriched than that for the thin disk, which is 
reflected in the Ba abundance while an enhanced SN II- 
like [a/Fe] is assumed for both. In the end, the adopted 
values of Z m are ([Fe/H], [Mg/H], [Ba/H]) = (-1.3, -0.9, 
-1.1), (-1.5, -1.1, -1.7) for thick and thin disks, respec- 
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tively. An initial gas of the thick disk formation is set 
to retain the same abundances as in Z m . The star for- 
mation rate (SFR) is assumed to be proportional to the 
gas fraction with a constant coefficient v for the dura- 
tion Agp. The higher value of z/=2 Gyr" 1 for Agp = 
1.5 Gyr is adopted for the thick disk, compared with 
v=0A Gyr -1 and Asf=12 Gyr for the thin disk. For 
the initial mass function (IMF), we assume a power-law 
mass spectrum with a slope x of —1.35 (a Salpeter's) to- 
gether with a mass ran ge (mi, m u )=(0.05 M©, 50 M©) 
(jTsuiimoto et al.lll997t) . 

In Figure 2, we show the model results for the evolu- 
tion of [Mg/Fe] and [Ba/Mg] against [Fe/H] in the thick 
disk. We see a good agreement with the observations. 
In particular, a successful reproduction of the [Mg/Fe] 
feature suggests that a new SN la DTD revealed by 
extragalaxy studies is compatible with the Milky Way 
case. In the [Ba/Mg]- [Fe/H] diagram, we find an ap- 
parent offset between in th e obs erved [Ba/Mg] va lues 
between IBensbv et all (j2005l ) and lVenn etafl (|2004D for 
bo th thick and thin disks. Therefore, the [Ba/Mg] data 
by IVenn et all (|2004h are shifted by +0.15 dex, which is 
equivalent to a mean difference between two data sets for 
-0.5<[Fe/H]<-0.2. 

We regard the thick disk as a first disk which is heated 
up by an ancient minor merger, that is subsequently fol- 
lowed by the gradual formation of a secondary disk, i.e., 
the thin disk. Such a first thick disk can also be formed 
through clump mergin g in an unstable primordial disk 
(|Bournaud et al.l I2007D . In these scenarios, we assume 
that star formation within the thin disk occurs after the 
termination of star formation in the thick disk, and thus 
the thin disk stars start forming from the thick disk's re- 
maining gas (corresponding to ~ 10% of the original gas) 
mixed with the gas accreted onto the disk. Accordingly, 
chemical abundances of the first stars in the thin disk are 
similar to those of the most metal-rich stars in the thick 
disk. As shown in Figure 2, the thick disk formation 
leaves the metal-rich gas as an end product of its chemi- 
cal evolution. As a result, the evolution of [Mg/Fe] starts 
from [Fe/H]~0 (Fig. 3). Then, [Fe/H] and [Mg/Fe] de- 
creases and increases, respectively, owing to dilution by 
metal-poor infalling gas. This reverse evolution comes to 
an end when the chemical enrichment by star formation 
exceeds the effect of gas dilution, and subsequently an 
usual evolutionary path appears. Our claim is that the 
remaining metal-rich gas after the thick disk formation 
results in the presence of no metal-poor thin disk stars. 
In the plot of [Ba/Mg], there are few observed data which 
coincide with the predicted early evolutionary path un- 
til [Fe/H]~-0.8. We note that it is not a contradictory 
result because the predicted fraction of stars riding this 
path is very small. 

4. CHEMICAL EVOLUTION OF BULGE 

We examine the chemical evolution of the Galactic 
bulge from a viewpoint that it is composed of two differ- 
ent populatio ns. Models are constructed based o n two- 
peaked MDF (IHill et alJl20TTt IBensbv et al.ll20Ll . Our 
grand picture of their formation is that (i) first, a metal- 
poor population (MPP) is formed with a relatively short 
timescale, and (ii) subsequently from its remaining gas 
plus a gas inflow from an inner disk, a metal-rich popu- 
lation (MRP) starts to form with a longer timescale. For 
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Fig. 4. — Predicted features of chemical evolution for the two- 
component bulge. The model results for metal-poor and metal-rich 
components are shown by blue and red curves, respectively. The 
calculated MDFs are convolved using the Gaussian with a disper- 
sion of 0.1 dex in [Fc/H] considering a measurement error expected 
in the data. The MDF summing each distribution up with a ratio 
of 0.5/0.5 is denoted by a blac k solid curve. The model distri- 
bution and the observed one by IBensbv et al.l l20TH l are normal- 
ize d to coinc i de wi th the total number of the sample stars used 
by IHill et al.l lIMUri .The observed dat a arc taken from Hill ct al. 
<|2011l filled circles) and IBensbv et alj l(20Tll_ open circles) for the 
MDF. IBensbv et all l|2011l , fille d circles) and Gonzalez et all (f20Tll . 
small crosses) for [Mg/Fe], and IBensbv et al.H2011l) for IBa/Mgl. 



the model of MPP, we assume a pre-enriched initial gas 
abundance, ([Fe/H], [Mg/H], [Ba/H])=(-1.3, -0.9, -1.1), 
so as to reproduce lack of stars for [Fe/H] <, -1.3 in the 
MDF and high Ba abundances of the most metal-poor 
stars. Abundances of an infalling gas are assumed to be 
the same. The enriched gas is considered to be an end 
result of chemical processing associated with the halo for- 
mation. The results calculated with (r; n , i>, Asf) = (0.3, 
4, 1) together with a Salpeter IMF are shown by blue 
curves in Figure 4. 

MRP is calculated with an initial condition given by 
a metal-rich gas abundance as an end result of a former 
MPP's formation under the setting of (r; n , v, Asf) = (1.5, 
3, 4). An inflow gas from an inner disk is assumed to 
be enriched up to ([Fe/H], [Mg/H], [Ba/H])=(-0.3, -0.3, 
-0.5). This subsolar metallicity is expected in the chemi- 
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cal evolution of the inner disk at an early epoch since at 
the inner disk an efficient chemical enrichment proceeded 
but the present-day gas abundance is not so metal- 
rich such as [Fe/H]~ + 0.2-0.3 inferred from Cepheids 
(jAndrievskv et aLl l2004) . This together with the consid- 
eration that sufficient time for the release of SN la and 
s-process products elapsed results in the above elemental 
abundances. In addition, in its modeling, we find that 
a top-heavy IMF is indispensable to make a metal-rich 
MDF as observed. Otherwise, a peak in the MDF results 
in being close to [Fc/H]~0 with an end metallicity un- 
reachable to [Fe/H]~+0.3. Here we assume x=l. 05. The 
results are shown by red curves. In the MDF (top panel), 
each contribution denoted by colored curves is summed 
up with an equal ratio according to an observed implica- 
tion, to construct the predicted MDF (black curve) that 
can be compared with the observed two-peaked one. We 
also see that the predicted [Mg/Fe] and [Ba/Mg] evolu- 
tions are consistent with the observed features. Note that 
a gas fraction at the end of star formation for each case 
of MPP and MRP is 9%, 8% respec tively in our calcu- 
lations. ICescutti fc Matteuccil (|201 lh first shows the Ba 
evolution in the bulge, calculating [Ba/Fe] together with 
the Mg/Fe evolution. Combination of these two evolu- 
tions broadly yields little evolution of [Ba/Mg] (~-0.4 — 
-0.5) at least for [Fe/H] £ 0. Its predicted level that is 
lower than the observed [Ba/Mg] on the whole is likely 
to be caused by a low Ba yield for the s-process in their 
models. 

Our new chemical evolution models have demonstrated 
that a top-heavy IMF is required to explain the ob- 
served MDF and dependences of [Mg/Fe] and [Ba/Mg] 
on [Fe/H] in the bulge stars. Some previous studies sup- 
port a top-heavy IMF in the bulge from a different an- 
gle. Ivan Dokkuml (|2008|) revealed that IMFs in early- 
type galaxies at 0.02 < z < 0.83 are signi ficantly natter 
than the present-day value of the Galaxy. iLarsonl (|1998f ) 
discussed a number of items of observational evidence 
that support top-hea vy IMFs in high-reds hift spheroidal 
galaxies. In addition. iManess et al.l(|2007|) claimed a top- 
heavy IMF in the Galactic center. We find an alternative 
model to explain the metal-rich part of MDF with the 
IMF unchanged. The device of this model is an assumed 
inflow which will chemically evolve to [Fe/H]~+0.5 dur- 
ing the bulge formation. However, its scenario implicitly 
presumes a top-heavy IMF in an inner disk outside the 
bulge. Therefore, we conclude that a top-heavy IMF 
is unavoidable for the chemical evolution of the inner 
Galaxy. 

In our scenario, metal-rich bulge stars are formed 
through an inflow from the inner thin disk. On the 
other hand, the bulge enriches the disk with large-scale 
winds to induce the pr oduction of metal-rich disk stars 
(jTsuiimoto et al.l [20101. In this way, we present a new 
view on Galactic chemical evolution that an interplay be- 
tween the bulge and the disk accelerates their chemical 
enrichment each other. 

5. DISCUSSION 

5.1. One- component bulge 

The complexity of the stellar population in the Galac- 
tic bulge has just entered the stage of debate on its poten- 
tiality. The one-component bulge is still a plausible view 
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Fig. 5. — Predicted features of chemical evolution including the 
MDF for the one-component bulge. The results with a natter IMF 
(:r=-1.05) are shown by solid curves while those with a Salpeter 
IMF (x=-l.35) are denote d by dashed curves. T he observed data 
for the MDF is taken from IFulbright et al.1 l2006h . 



in ter ms of the MDF (IFulbright et al.ll2006tlZoccali et all 
I2008H or stellar kinematics (|Howard et aLl l2009t). ThusT 
it is worthwhile to refer to the chemical evolution based 
on a single-pea ked MDF derived from red giants in the 
Baade window (|Fulbright et al.l [2006f) . The reason dif- 
ferent typ es of stars yields di fferent MDFs remains un- 
explained (Be nsbv et alJl20lTl for a detailed discussion). 
First, we show that the model with an enhanced star 
formation is not sufficient to reproduce the chemical fea- 
ture of bulge stars. The dashed curve in the top panel of 
Figure 5 represents a resultant MDF calculated by the 
model with a rapid collapse (rj n =0.3 Gyr), a high SFR 
(u=2 Gyr -1 ), and Asf=2 Gyr together with a Salpeter 
IMF. These setting results in a final gas fraction of 10%. 
For the metallicity of an initial gas as well as of an in- 
falling gas, we assume a very low-metal content, ([Fe/H], 
[Mg/H], [Ba/H]) = (-3.0, -2.6, -3.2). The MDF thus ob- 
tained is entirely skewed to a low metallicity. In addition, 
the predicted [Mg/Fe] curve in the middle panel is lower 
than the observed data in a metal-rich regime. 

On the other hand, the predicted [Ba/Mg] in the bot- 
tom panel exhibits a sharper rise from a much lower 
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metallicity than is expected from the observation. These 
inconsistencies are resolved by the models with a flatter 
IMF (x—1.05), as shown by solid curves, which give a 
good agreement with the observed MDF as well as the ob- 
served correlations of [Mg/Fe] and [Ba/Mg] with [Fe/H]. 
This flat IMF is fairly consistent with those claimed by 
sever al authors, such as x = 1.1 (|Matteucci fe Brocatol 
119901 ) and x = 0.95 (Ballero et al. 2007a,b). The neces- 
sity of a top-heavy IMF in the bulge can be also con- 
cluded from the following argument. An anticipated 
rapid collapse realizes the condition approximated by 
a closed-box model. In this approximation, the mean 
metallicity of stars will become the heav y- element yiel d 
when the remaining gas approaches zero (Tinslcy 1980). 
In other words, the location of a peak in the MDF is 
fundamentally determined by the IMF (not by the SF 
efficiency) . 

In the end, our study has demonstrated that a top- 
heavy IMF is the key factor to efficiently drive chem- 
ical enrichment in the bulge, whether one-component 
or two-components. Recently revealed two-component 
bulge should be validated by a bigger database of a large 
sample of bulge stars to confirm the dip around [Fe/H]~0 
in the MDF. 

5.2. Halo vs. short- delayed SNe la 

The well-known fact for the Galactic halo is basically 
no indication of SNe la for the elemental abundance of 
halo stars, which exhibits a plateau of [a/Fe] ratio over 
a whole metallicity range. It turns out that in the CDM 
Universe halo stars were rapidly formed in the Galac- 
tic building blocks, likely in dwarf galaxies, with its ter- 
mination probably due to huge energy released by nu- 
merous SNe II before the major occurrence of of short- 
delayed SNe la. Thus, individual building blocks must 
be formed with a very short timescale (~ 10 8 yr), while 
an assembly of them finally make s the stellar halo whic h 
exhibits an age span of a few Gyr (|De Angeli et al.ll2005h . 
Observational studies reveal that the chemical abun- 
dances of the Galactic stellar halo are significantly differ- 
ent from those of the present-day dwarf galaxies around 
the Galaxy in the sense that halo stars have higher |al- 
pha/Fe] (e.g.. IShetrone et al.l 120011: iTolstov et ail 120031: 

I Venn et al.N2004ft . Recent theoretical calculations based 
on the A CDM model show that the majority of the 
Galactic halo stars are formed from a few relative mas- 
sive dwarfs in which star formation is rather rapid and 
thus the stars are chemically enriched primarily by SNe 

II to have an enhanced [alpha/Fe] ratio (|Robertson et al.l 
2005). These numerical results should be reanalyzed by 
introducing the new DTD for SNe la. 

It would be interesting to check if the s-process el- 
ements from AGB stars starting to release with a 
timescale of a few 10 8 yr are imprinted in the abundances 
of halo stars. In the case where the s-process operated 
in the stellar halo, the [s-process/r-process] ratio such 
as [Ba/Eu] or [La/Eu] switches its constant value to an 
increasing feature from a certain metallicity owing to a 
gradual s-process contribution superimposed on the r- 
process material in the ISM. Observationally, the metal- 
licity indicating this onset among the Galactic halo stars 
is difficult to detect si nce it is hidden by a large scatter in 
the abundance ratios. iGilrov et al.l (|1988f ) find the onset 
of s-processing at [Fe/H]~-2. On the other hand, recent 



study claims no indication of the s -process at least until 
[Fe/H]~-1.4 ([Roederer et al.llMfl . For [Fe/H]>-1.4, ei- 
ther the presence or the absence of s-processing is very 
unclear because [La/Eu] continues to make a plateau 
while [Pb/Eu] shows an upward trend with an increasing 
[Fe/H]. This issue points to the need for further studies. 

If the signature of s-processing will be detected among 
halo stars, how do we connect the new DTD to the chem- 
istry of halo stars? A deficiency of binary stars in the 
Galactic globu lar clusters (GCs) compared to the field is 
reported (e.s., iPrvor et al.lll989l; [Cote et al.lll996t) . Sup- 
pose that this observed fact results from the lack of bina- 
ries with short separations in GCs, some of which will be 
the progenitors of short-delayed SNe la in either a DD 
scenario or a SD scenario. Then, the GC stars are ex- 
pected to be hardly enriched by SNe la. These GCs are 
suggested to be originally form ed within galactic build- 
ing blocks (Se arle fe Zindll978l ). and halo stars can be re- 
garde d as the stripped stars of the blocks (Bck ki fe Chibal 
2001 s ). Therefore, we can hypothesize that in individual 
blocks both GCs and halo stars shared giant molecu- 
lar clouds, where a lower fraction of close binary stars 
could be formed due to their high densities, as their 
birth place. In addition, recent observed result based on 
chemically unusual red giants implies that a significant 
fracti on of halos stars ar e originated from disintegrated 
GCs (|Martell et al.ll2011| ). Accordingly, halo stars could 
avoid the contamination by SNe la in their chemical 
history. The compatibility of the presence of s-process 
among halo stars with the short-delayed DTD can be 
also understood consistently if s-proce ss elements are 
produced in fast-rotating m assive stars (jpignatar i et al.l 
[20081: iChianoini et al.ll20Ll . 

6. CONCLUSIONS 

The new model for the Galactic bulge with two 
episodes of star formation is presented. We show that it 
explains the two-peaked MDF recently deduced from red 
clump stars, that can avoid a contamination of disk stars 
more reliably than red giants, as well as from microlensed 
dwarf stars. We stress that this successful MDF repro- 
duction is accompanied by a coincidence of the observed 
and predicted evolutions of [Mg/Fe] and [Ba/Mg]. In 
particular, we first examine the [Ba/Mg] evolution in the 
bulge, adopting a new Ba yield, which results in the com- 
patibility with an observed level of Ba abundances for 
both disks and the bulge. Our result will give feedback 
to the s-process nucleosynthesis calculations. We find 
that an early [Ba/Mg] trend supports a pre-enriched gas 
for the bulge formation, which is consistent with lack of 
metal-poor stars in the observed MDF with two peaks. 

The new DTD for SNe la is found to be compatible 
with the chemical evolution of the Galaxy, leaving room 
for further discussion on chemical abundances of halo 
stars. This short-delayed SNe la virtually enhance the 
Fe production with a short timescale together with SNe 
II. Our study, however, confirms that a top-heavy IMF 
is still required to explain the chemical feature of the 
Galactic bulge. Many works thus far have been devoted 
to the investigation of the effect of various SN la DTDs 
on chemical evolution in the sol ar vicinity, focusing on 
the observed [a/Fe] break (e.g., iMatteucci et al.l 120091 : 
iKobavashi fe Nomotol [20091 . On this matter, we claim 
that the discussion on the DTD in this context should 
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be assessed by comparing the modeled chemical feature 
of the thick disk with the corresponding observed one. 
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